First order phase transition in the frustrated triangular antiferromagnet CsNiCls. 



NO 
O 
O 
(N 

G 

(N 



o3 



C 

O 

o 



> 

ON 

o 

NO 

o 

NO 

o 

-I— > 
03 



T3 

o 
o 



X 



G. Quirion, X. Han, M.L. Plumer 1 and M. Poirier 2 

Department of Physics and Physical Oceanography, 
Memorial University, St. John's, Newfoundland, Canada, A1B 3X7 
2 Regroupement Quebecois sur les Materiaux de Pointe, Departement de Physique, 
Universtte de Sherbrooke, Sherbrooke, Quebec, Canada, 31K 2R1 
(Dated: February 6, 2008) 

By means of high-resolution ultrasonic velocity measurements, as a function of temperature and 
magnetic field, the nature of the different low temperatures magnetic phase transitions observed for 
the quasi-one-dimensional compound CsNiCl3 is established. Special attention has been devoted to 
the field-induced 120° phase transition above the multicritical point in the H — T phase diagram 
where the elastic constant C44 reveals a step-like variation and hysteresis effects. These results 
represent the first experimental evidence that the 120° phase transition is weakly first order and 
contradict the popular notion of new universality classes for chiral systems. 



Despite increasing interest in recent years in geometri- 
cally frustrated systems , there still exists considerable 
controversy over the nature of the phase transition in 
the prototypical magnetic system characterized by near- 
neighbor antiferromagnetic exchange interactions on a 
stacked triangular lattice 0, 0, Q . The prediction made 
twenty years ago that helical degeneracy associated with 
the 120° magnetic order leads to new Heisenberg and 
XY chiral universality classes [f| found support over the 
following ten years or so from numerous renormalization- 
group studies, Monte-Carlo simulations and experimen- 
tal data 0, IH El 0- Results in favor of this scenario 
continue to appear j^i fiol . An alternative proposal 
of a very weak fluctuation-induced first order transition 
made soon after the original suggestion nas a l so been 
strengthened by further theoretical studies and numerical 
simulations |l3lll4ll5Ulfill7lll8| . Experimental data on 
both rare-earth helimagnets as well as ABX3 compounds, 
such as CsMnBr3, have been used extensively to support 
each scenario. Evidence for a weak first-order transition 
was found in the thermal expansion data on the helimag- 
net Ho long before this controversy surfaced The 
results presented in this Letter reveal the first experi- 
mental evidence that the 120° XY transition is weakly 
first order. This is achieved also through magnetoelastic 
coupling effects, via ultrasonic sound velocity measure- 
ments in the field-induced 120° phase of CSNICI3. 

CsNiCbj is one of the more widely investigated of 
the large class of quasi-one-dimensional hexagonal ABX3 
materials with strong antiferromagnetic c-axis exchange 
0. Along with its sister compounds CsMnI.3, CsNiBr3 
and RbNiBr3, it has weak c-axis anisotropy giving rise 
to a Linear (L) Ising-like ordered state in zero field at 
Tjvj = 4.8 K, with an additional in-plane ordering at 
Tjv 2 = 4.4 K, resulting in an elliptical (E) polariza- 
tion of the spin density (S) at low temperatures [20j| . 
These phases are characterized by a period-2 modula- 
tion along the c-axis and period-3 in the basal plane. A 
magnetic field applied along the c-axis of only 2.3 T is 
sufficient to induce a spin-flop phase where S now lies in 



the basal plane and forms the familiar 120° spin structure 
of the frustrated triangular antiferromagnet. These three 
phases meet at an unusual type of multicritical point |2l| 
(at T m = 4.53 K, H m = 2.3 T). The transition to the 
Ising-like state involves a phase-factor degeneracy asso- 
ciated with the triangular geometry and is been predicted 
to belong to the XY universality class. The nature of even 
this transition has a history of controversy, with the most 
recent addition being an analysis of neutron diffraction 
data on a similar transition in CsCoBr3 suggesting tri- 
critical behavior |22j . The transition at Tn 2 is generally 
accepted to also be of XY universality. At the multi- 
critical point, the axial anisotropy is exactly canceled by 
the applied field and the system is effectively isotropic 
(Heisenberg-like) . At higher field strengths, there is XY 
symmetry along with chiral degeneracy. Such systems 
provide a convenient platform for the investigation of 
a line of phase transitions (to the paramagnetic state) 
which can be sampled by changing the field strength. In 
contrast with present data, previous experimental inves- 
tigations of this transition boundary support the notion 
of n — 2 chiral universality 0, 0| • 

Magnetoelastic coupling has been repeatedly demon- 
strated to be a useful mechanism to reveal the nature of 
the magnetic ordering in CsNiCls j^, [2^, ■ Landau 
theory can be used to show how elastic constants scale 
with the various order parameters and also to yield mean- 
field predictions of anomalies at the transition boundaries 
|28j . In the present work, results and analysis are pre- 
sented of high-resolution measurements of various elastic 
constants as a function of both temperature and mag- 
netic field, with a focus on the paramagnetic-to-120° 
spin phase boundary. Step-like discontinuities are found 
where Landau theory predicts none. Attempts at curve 
fitting to extract a critical exponent [3 show that this 
leads to values which are field dependent. The strongest 
evidence that the paramagnetic to spin-flop phase bound- 
ary is weakly first order is found in the hysteretic behav- 
ior of C44 as a function of temperature and field. 

Contrary to previous investigations [2^, I27I l2flj , where 
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ultrasonic techniques were mainly used in order to deter- 
mine phase diagrams, the emphasis of the present work 
is on the measurement of critical phenomena in CSN1CI3. 
This was realized using a high-resolution pulsed ultra- 
sonic interferometer to measure the temperature and 
magnetic field dependence of different acoustic modes 
propagating along or perpendicular to the hexagonal c- 
axis. Measurements were carried out on a single crystal 
specimen of 8.9 mm in length along the c-axis and ap- 
proximately 2.5 mm along the perpendicular directions. 
The acoustic modes at 30 MHz were generated using lon- 
gitudinal and transverse lithium niobate transducers. 

The analysis of the temperature or field dependence 
of elastic properties provides a convenient way to study 
magnetic critical behavior in many systems. For hexag- 
onal CSN1CI3, a Landau type approach has been used to 
determine the relationship between the variations in the 
elastic constants C33 and Cn and the various order pa- 
rameters [28J. This coupling occurs due to magnetoelas- 
tic contributions to the free energy of the form ~ g e i S 2 , 
where is an element of the strain tensor and S is the 
order parameter. A similar approach can be used to in- 
clude coupling terms which account for the application of 
a magnetic field along the c-axis. One result of this new 
model 30] indicates that the relative variation AC33 / C33 
can be generalized as 
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where A and 7 are constants specific to the magnetic 
transition of interest, while the temperature and field 
dependencies are directly associated with those of the 
order parameter S. According to QJ, C33 is expected 
to show a discontinuity even in the case of a continuous 
phase transition. Thus, results on C33 cannot be used to 
discriminate between a continuous and weakly first or- 
der transition. This type of behavior is to be expected 
whenever a linear-quadratic coupling, between the strain 
and the order parameter, is allowed by symmetry |28| . 
In order to clearly determine the character of the tran- 
sition of the 120° phase, other ultrasonic modes need to 
be used, in particularly those that depend exclusively 
on quadratic-quadratic coupling terms (e 2 S 2 ). The al- 
lowed coupling terms, compatible with hexagonal sym- 
metry (not included in our previous analysis [2^]) are 
simply 



F c (S,e t ) = ^(e 2 + e 2 )S 2 + 9 -fe 2 S 2 z 
+ 9 -^{el + el)Sl + 9 -^elSl 



(2) 



where the notation of Ref. [28] has been used. As these 
coupling terms are quadratic in strain, the variation in 
the elastic constants can be written as 
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In the case of an hexagonal structure, C44 and Cqq can be 
obtained by measuring the velocity of transverse waves 
propagating along and perpendicular to the c-axis, re- 
spectively. 

Figure ^ shows the most significant results of the tem- 
perature dependence of AC66/C66 with the magnetic 
field applied along the c-axis. At H = T, the onset of 
the L - E phase transition is clearly visible (Tjv 2 = 4.33 K) 
while the variation at is barely noticeable. The re- 
sults show two distinct behaviors depending on whether 
the value of the field is lower or higher than H m . In 
the elliptical phase, Cqq softens as the temperature de- 
creases while the opposite trend is observed in the 120° 
spin phase above H m . The observed temperature depen- 
dencies of AC66/C66 are perfectly consistent with the 
Landau predictions and these data can be used to es- 
timate the temperature dependence of the order parame- 
ter S±. The results of this analysis, obtained at different 
fields, are presented in Fig.|21as a function of the reduced 
temperature r = 1 — T/Tn 2 on a log-log plot. All curves 
show a well defined power law behavior over a minimum 
of two decades in r. Clearly, for H < 1.5 T a unique 
scaling is observed, confirming that (3e = 0.35 ± 0.02 is 
field independent in the elliptical phase. As the value 
of the field approach H m , the value of the critical expo- 
nent (3 suddenly decreases and then gradually increases 
at higher fields. The insert in Fig. illustrates in detail 
how the value of f3 evolves as a function of a magnetic 
field for CSN1CI3. Very close to the multicrital point, j3 
reaches a minimum of (3 = 0.25 ± 0.02 which corresponds 
to the predicted value for n — 2 chirallity 5] but is also 
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FIG. 1: Relative variation of the elastic constant Ce6 as a 
function of temperature. The broken and continuous lines 
represent results obtained below and above H m = 2.3 T, re- 
spectively. 
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FIG. 2: Order parameter as a function of the reduced tem- 
perature t = 1 — T/Tjv 2 calculated using J1J and the results 
presented in Fig. . The data obtained at different field 
are represented by symbols while lines represent fitted data 
at small r using a simple power law. The inset shows the 
field dependence of the critical exponent j3 obtained from the 
power law fits. 



close to that of tricriticality, /? = 1/4. At higher fields, 
P increases significantly. As noted previously 0, Il7| . 
variation in the value of effective critical exponents as a 
function of an irrelevant parameter (in this case, H) may 
indicate a weak fluctuation-induced first order transition. 

Close inspection of the upper curves presented in Fig.Q] 
show an unexpected dip right at the para-120° phase 
boundary. This very small anomaly persists at all fields 
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FIG. 3: Relative variation of the elastic constant C44 as a 
function of temperature. The broken and continuous lines 
represent results obtained below and above H m = 2.3 T, re- 
spectively. 



above H m and could be interpreted as the effect of a 
weakly first order transition. Motivated by these re- 
sults, additional measurements were made using trans- 
verse waves propagating along the c-axis and polarized 
in the basal plane to obtain C44. This series of results 
was obtained with an exceptionally high resolution of 
0.1 ppm. The data are presented in Fig. [3] as a func- 
tion of temperature. At H = 0, two distinct features 
are noticeable and correspond to the onset of the phase 
transitions at T^ t and Tjy 2 . As the field is increased, 
both anomalies merge at H m = 2.3 T. Above the multi- 
crital point (see continuous lines), the observed temper- 
ature dependence changes within the 120° phase as the 
field increases. Moreover, no power law relationship, as 
predicted by the Landau Model J3J), could be identified. 
More significantly, a step like variation is noticeable at 
the critical temperature. These observations taken all 
together cannot be reconciliated with a continuous phase 
transition. 

As a test of the first order character of the para- 120° 
phase transition, the possibility of thermal hysteresis was 
investigated. A collection of thermal cycles, realized us- 
ing a cooling/heating rate of 0.1 K/min, is presented in 
Fig. 0] The results obtained for the 120° phase are com- 
pared to the data collected at H = T. At zero field, 
where all experimental evidence presented in this paper 
clearly indicate that the phase transition is continuous, 
no significant hysteresis is observed. However, the transi- 
tion to the 120° spin phase show a difference between the 
data collected during the heating and cooling processes. 
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FIG. 4: Thermal cycle analysis realized at different fields 
using the relative variation of the elastic constant C44. The 
data collected during the cooling process are all represented 
in red color for clarity. Curve (a) corresponds to the zero 
field transition at T/v 2 . All other curves are associated with 
the transition to the 120° phase. 
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FIG. 5: Relative variation of the elastic constant C44 as a 
function of H 2 measured at T = 5 K. The inset shows the 
relative variation of C44 after subtracting the magnetoelastic 
contribution observed in the paramagnetic phase. 



Those differences are small but are systematically ob- 
served at all field values. The hysteresis is maximum just 
below the critical temperature and persist over a temper- 
ature range of about 0.5 K. An additional confirmation 
of first-order character is obtained from the field depen- 
dence of C44. The data collected at T — 5.0 K, presented 
in Fig. as a function of H 2 , clearly show the quadratic 
field dependence of C44 in the paramagnetic phase. This 
field dependence is associated with magnetostriction ef- 
fects previously observed in CSN1CI3 31]. The field de- 
pendence of C44 near the 120° phase transition shown in 
the insert of Fig. [S] has been isolated by subtracting this 
magnetostriction. It is clear that the step like variation 
of C44, along with the observed hysteresis, at the 120° 
phase boundary cannot be accounted for in the context 
of a continuous phase transition. 

The data and analysis presented in this work serve 
to fill a long-standing gap regarding experimental sup- 
port for the growing body of theoretical and numerical 
work that the phase transition in the prototypical geo- 
metrical frustrated system, the stacked triangular anti- 
ferromagnet, is fluctuation-induced first order in nature. 
CSN1CI3 provides a convenient field-temperature phase 
diagram for this purpose with a phase boundary line to 
the 120° spin structure which may be sampled at vari- 
ous points. Sound velocity measurements have proven to 
be an accurate tool for obtaining high resolution data on 
the various order parameters. The possibility to extract 
effective (field dependent) critical exponents, the weak 
nature of the discontinuities and small hysteresis, taken 
together, provide evidence for the weakness of the first- 
order character of this transition. It is thus not surprising 
that conventional renormalization-group techniques and 



Monte-Carlo simulations had previously been supportive 
of the notion of new chiral universality classes associated 
with such frustrated systems. 
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